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ABSTRACT 

We present five new absolutely calibrated continuous stellar spectra constructed as far as possible from 
spectral fragments observed from the ground, the Kuiper Airborne Observatory {KAO), and the IRAS Low 
Resolution Spectrometer. These stars — a Boo, y Dra, a Cet, y Cru, and p UMa — augment our six, 
published, absolutely calibrated spectra of K and early-M giants. All spectra have a common calibration 
pedigree. A revised composite for a Boo has been constructed from higher quality spectral fragments than 
our previously published one. The spectrum of y Dra was created in direct response to the needs of 
instruments aboard the Infrared Space Observatory (ISO)\ this star’s location near the north ecliptic pole 
renders it highly visible throughout the mission. We compare all our low-resolution composite spectra with 
Kurucz model atmospheres and find good agreement in shape, with the obvious exception of the SiO 
fundamental, still lacking in current grids of model atmospheres. The CO fundamental seems slightly too 
deep in these models, but this could reflect our use of generic models with solar metal abundances rather 
than models specific to the metallicities of the individual stars. Angular diameters derived from these 
spectra and models are in excellent agreement with the best observed diameters. The ratio of our adopted 
Sirius and Vega models is vindicated by spectral observations. We compare IRAS fluxes predicted from our 
cool stellar spectra with those observed and conclude that, at 12 and 25 ju,m, flux densities measured by 
IRAS should be revised downwards by about 4.1% and 5.7%, respectively, for consistency with our absolute 
calibration. We have provided extrapolated continuum versions of these spectra to 300 pm, in direct support 
of ISO (PHT and LWS instruments). These spectra are consistent with IRAS flux densities at 60 and 100 
pm. © 1996 American Astronomical Society. 


1. INTRODUCTION 

The previous papers of this series have presented a con- 
sistent effort to provide^ absolutely calibrated broadband and 
narrowband infrared photometry and spectra based upon a 


carefully selected, infrared-customized pair of stellar models 
for Vega and Sirius, created by Kurucz, and absolutely cali- 
brated by Cohen et al. (1992a, hereafter referred to as Paper 
I). These hot stellar models have been employed as reference 
spectra to provide continuous (i.e., uninterrupted in wave- 
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length), observed, absolutely calibrated spectra (“compos- 
ites”) of cool stars by the methods detailed by Cohen et al. 
(1992, hereafter referred to as Paper II) and Cohen et al. 
(1995, hereafter referred to as Paper IV). This approach has 
yielded six infrared-bright secondary stellar standards with 
calibration pedigrees directly traceable to our primary radio- 
metric standard, namely, a CMa. Molecular bands of CO 
and SiO dominate these observed infrared spectra of cool 
giants and supergiants (Cohen et al. 1992b, hereafter referred 
to as Paper III). 

In the present paper we update the spectrum of a Boo to 
take account of newer and higher quality spectral fragments, 
and we augment the six composites by spectra of four addi- 
tional stars assembled by the identical procedures to those 
described in Papers II and IV (Sec. 2). Note that two of these 
new composites do not start at 1.2 /x m because there has 
never been any further development of the unique Strecker 
et al. (1979, hereafter referred to as SEW) archive of air- 
borne 1.2-5. 5 /xm spectra. The composite of /x UMa starts at 
about 3.1 (x m, and that for -y Cru at about 3.9 /xm. However, 
these can still provide valuable calibrators (cf. Cohen et al. 
1996a, hereafter referred to as Paper VI which offers 3.0—29 
/xm coverage). 

In a future paper (Cohen et al. 1996b, hereafter referred to 
as Paper VIII), we will describe how we have converted this 
limited set of stars with abundant spectral observations into a 
network of nearly 300 calibrators across the sky by using the 
“template assumption,” namely, that the dereddened infra- 
red spectrum of any observed K-M giant accurately repre- 
sents the intrinsic spectrum of any other giant with the iden- 
tical spectral type as the composite, and so serves as a 
‘ ‘template. ’ ’ Given our intended dependence upon these few 
observed continuous spectra (now covering the range K0 
—M3. 4 III), it is appropriate to make direct comparisons 
with a grid of model atmospheres. We, therefore, compare 
Kurucz (1991, 1993) models for our stars with composite 
spectra, and ratio these calibrated spectra to the models (Sec. 
3). 

We have several objectives in these comparisons: to test 
the goodness of fit (shape) of theory and observations (these 
cool stellar models are totally independent of the data and 
procedures used to assemble the composites); if well 
matched in shape, to determine angular diameters from the 
normalization of models that match the composites; to com- 
pare these empirical angular diameters with those observed 
by Michelson (intensity) interferometry or lunar occultation; 
and to attempt the extrapolation of these observed spectra to 
the far-infrared, using IRAS measurements to attempt to vali- 
date these extrapolations. Despite the conspicuous omission 
of the SiO fundamental (whose presence locally dominates 
these spectra between about 7.5 and 12 /xm), we find that the 
Kurucz models are capable of replicating these complete 
spectral shapes to within a few percent. Our derived angular 
diameters are in good accord with the best, observationally 
determined, diameters for these stars (Sec. 4). 

We have sought direct observational evidence on the ratio 
of our adopted spectra of Sirius and Vega, both photometri- 
cally and through 10 /xm spectroscopy, and find support for 
this ratio (Sec. 5). 


Table 1. New composite spectra described in this paper.* 


Star 

Spectral type 

Date of assembly 

a Boo 

K1 III 

December 21, 1995 

7 Dra 

K5 III 

December 22, 1995 

a Cet 

Ml. 5 III 

March 26, 1996 

7 Cru 

M3. 4 III 

March 13, 1996 

fi UMa 

M0 III 

April 2,1996 


*Table 1 can be found in the AAS CD-ROM Series, Vol. 7, 1996. 


In Sec. 6, we return to the issue of the mid-infrared be- 
havior of real giants, as opposed to model atmospheres, and 
present empirical evidence in favor of the locally Rayleigh- 
Jeans character of K and M giants beyond 10 /xm. 

Finally, we again probe the absolute calibration of IRAS 
(cf. Paper I) by synthesizing IRAS flux densities for all our 
composites in the 12 and 25 /xm bands (Sec. 7) and, from 
theoretical continuum extrapolations of their spectra using a 
grid computed by one of us (D.F.C) (Sec. 8), in the 60 and 
100 /xm bands (Sec. 9). We find good agreement between 
these extrapolated models and observed IRAS long- 
wavelength flux densities, encouraging us to provide stellar 
calibration spectra in support of ISO’s absolute calibration. 

2. NEW COMPOSITE SPECTRA 

Our methodology for creating a complete and continuous 
composite spectrum from 1.2 to 35 /xm follows that de- 
scribed in Papers II and IV. Except for that of y Cru, these 
new complete spectra were made possible through a fruitful 
collaboration with the UKIRT Service Observing program 
(cf. Cohen & Davies 1995, hereafter referred to as Paper V) 
that has provided high quality 10 and 20 /xm spectral frag- 
ments. Table 1 summarizes the new composites. Full details 
of the process of assembly for each composite appear in the 
AAS CD-ROM Series, Vol. 7, 1996 associated with this pa- 
per, including the independent photometry used to support 
the normalization of spectral fragments. To illustrate the cali- 
bration pedigree that accompanies each new composite spec- 
trum, we present only one table here, corresponding to a 
Boo. Figure 1 presents the spectra themselves in the form of 
log X 4 jF x vs log X plots. We show the three composites with 
fullest wavelength coverage before those with more re- 
stricted ranges. 

All the new spectral fragments incorporated into this pa- 
per were secured with two or more of the following instru- 
ments: the NAS A- Ames “SIRAS” (Short Wavelength Infra- 
red Array Spectrometer: Wooden 1989), “FOGS” (Faint 
Object Grating Spectrometer: Witteborn & Bregman 1984), 
or “HIFOGS” (High-efficiency Infrared Faint Object Grat- 
ing Spectrometer: Witteborn et al. 1995) on the KAO or with 
the NASA 1.5 m Mt Lemmon telescope; the CGS3 spec- 
trometer on the 3.8 m UKIRT; or the IRAS Low Resolution 
Spectrometer-(hereafter referred to as LRS). The Ames spec- 
trometers are doubly sampled by using two different grating 
settings spaced an integral number plus one half detectors 
apart to achieve Nyquist sampling and to provide coverage 
of occasional dead detectors. CGS3 is always either doubly 
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Fig. 1 . The five cool stellar composites in the form of logarithmic plots of 
absolute \ 4 f\ against X. In each spectrum, the mean ±ler error bounds are 
indicated by the pair of solid lines. Short dashed line is the best-fitting 
Engelke approximation, (a) a Boo, (b) y Dra, (c) a Cet, (d) y Cru, (e) 
UMa. In all these plots, we represent the photometry points used to build the* 
composite by filled squares. Other photometry, unused by our procedures, 
can be used to validate the spectra, and such points are shown by open 
squares (ground-based data) or by open triangles ( IRAS data). 


0 .5 1 1.5 

(c) lo gto M/^m) 



or triply sampled, likewise to achieve the full theoretical 
resolution. For the CGS3 fragments (cf. Paper V) we have 
continued to emphasize close matching of airmasses between 
the cool target stars and the reference object (preferably 
Sirius, or one of our already established cool composites). 
From the KAO we could not always exert such control over 
stellar airmasses at the time of observation due to the relative 
brevity and constraints of the actual observing time at alti- 


tude. Therefore, we remove the residual effects of the terres- 
trial atmosphere (particularly important from 5 to 8 /urn) post 
facto from the spectra, through use of precomputed standard 
atmospheric transmission curves relevant to different 
amounts of precipitable water vapor (currently between 6 
and 12 yam). Actual water vapor was measured on each flight 
with a radiometer. These transmissions were calculated using 
Lord’s (1992) “ATRAN” software tool. Any corrections to 
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Table 2. Engelke functions used in, or relevant to, composite spectra.* 


Star 

Effective Temperature 

(K) 

Angular diameter 
milliarcsec 

a Boo 

4362 

21.12 

7 Dra 

3986 

10.17 

a Cet 

3745 

12.77 

7 Cru 

3626 

26.14 

fx UMa 

3735 

8.32 


*This table can be found in the AAS CD-ROM Series, Vol. 7, 1996. 


measured spectral points were ~3%-5% in molecular lines. 

As in Paper IV, we have substituted the best-fitting En- 
gelke (1992) function for noisy LRS, 16-24 /xm CGS3, 
and/or 16-35 /tm KAO data to provide meaningful extrapo- 
lations of our new composites out to 35 /zm. This analytic 
approximation to real stellar continua is derived assuming 
that the dominant source of infrared continuum opacity in 
these cool giant stars is H free-free (cf. Engelke 1992). 
There is only one essential input to the shape of an Engelke 
function, namely the effective temperature, which we fixed 
by reference to the literature (see Tables 2, 3, and the AAS 
CD-ROM Series). The normalizing scale factors for these 
functions, namely their angular subtenses, are determined by 
the fit of the Engelke function shape to our composite spec- 
tra. 

We document the full calibrational information for these 
new composites (see Table 4 and the CD-ROM), namely, 
date of assembly (i.e., the version of this spectrum); details 
of the photometry used to calibrate the composite radiometri- 
cally including the FWHM of the relevant passbands and 
their effective wavelengths for the Vega spectrum and for the 
star in question, and the monochromatic flux densities; all 
archival spectral fragments, with total spectral ranges, ranges 
actually utilized, and average resolving power over that spec- 
tral range (expressed as the resolving power, A/AX). We had 
access to adequate photometry to assemble the composites 
for a Boo, y Dra, fi UMa, and y Cru by precisely the same 
technique used for that of a Tau (Paper II) and the 1.2-35 
/urn. spectra detailed in Paper IV. For a Cet we lacked any 
traceable photometry in characterized passbands in the 10 
Atm window and were, therefore, obliged to assemble this 
star’s composite using the “end-to-end” method used in Pa- 
per VI for the composite spectrum of a 1 Cen. 

The star /x UMa is now designated as a new suspected 
variable (NSV 4829) by Kholopov et al. (1992), although 
details are minimal. Hoffleit & Warren (1991) indicate it to 
be a spectroscopic binary, possibly an eclipsing variable, 
with visual amplitude about 0.3 mag. Examination of the 
infrared literature (Gezari et al. 1993) does not suggest in- 
frared variation above a few hundredths of magnitude, after 
allowance for the diverse zero points of the references cited 
therein, although Sinton & Tittemore (1984) indicate quali- 
tatively that they saw some variation at M. Consequently, we 
have retained this star as a secondary standard because of its 
popularity as an-infrared calibrator, due as much to its north- 
erly location as to its mid-infrared brightness. 

In addition to the photometry used for the actual assembly 
of composites (the filled squares in Fig. 1), we sometimes 


have other photometry in well-characterized passbands that 
can be used to validate the spectra, precisely because it was 
unused by our procedures. Any such points are distin- 
guished, in Fig. 1 , by open squares (ground-based data) or by 
open triangles ( IRAS data, similarly not used in composite 
assembly). All photometry points are plotted at their isopho- 
tal wavelengths for these stellar spectra. These validating 
data points are as follows. For a Boo we have: narrowband 
Jn (Selby et al. 1988; Hammersley 1995, described in Paper 
VIII); UKIRT broad A- band (10 /zm) and Q -band (20 /zm) 
points (too wide to be integrated over the CGS3 10/20 /zm 
fragments but valid post factor; and the IRAS FSS flux den- 
sities at 12 and 25 /zm. For y Dra, these are narrowband Jn 
(Selby et al. 1988; Hammersley 1995) and broadband Ten- 
erife J (Hammersley 1995); a UKIRT Q-band point; and 
IRAS FSS points. For the other three stars, we have only 
IRAS flux densities as external checks of the composites. In 
all five cases, these additional photometry values confirm the 
accuracy of the new composites at the 1 cr level or better. 

We prefer to provide pristine data whenever possible, 
rather than to regrid each composite to some equally spaced 
or common wavelength scale. Each composite (on the CD- 
ROM), therefore, has a different set of wavelengths. We 
tabulate: wavelength; monochromatic flux density (F k in 
units of W cm _2 /zm _1 ); total uncertainty (also in units of 
W cm' 2 /im -1 ) associated with this value of F k ; local bias; 
and global bias. For most applications, “total uncertainty” is 
the error term most appropriate to use. It is the standard 
deviation of the spectral irradiance and includes the local and 
global biases. Local and global biases are tabulated as per- 
centages of the irradiance; their corresponding absolute 
quantities are already included in the total uncertainty. The 
global bias does not contribute error to flux ratios or color 
measurements, and may be removed (in the root-sum-square 
sense) from the total error. 

3. COMPARISON WITH KURUCZ MODEL ATMOSPHERES 

De Jager & Nieuwehuisen (1987) tabulate statistical rela- 
tions between effective temperature and bolometric luminos- 
ity for all spectral types, assigning algebraic variables, that 
can be interpolated, to the qualitative descriptors of “spectral 
class” and “luminosity class.” Following the approach and 
work of de Jager & Nieuwehuisen, Cohen (unpublished) has 
constructed similar tables that correspond to stellar mass and 
radius, using values drawn from the literature, whence a de- 
rivative relation between gravity and spectral type was ob- 
tained. Adopting solar metallicities, appropriate effective 
temperatures, and these gravities (corresponding solely to the 
spectral types of our stars), we have extracted theoretical 
spectra from the Kurucz grid. To provide a generic model for 
each selected set of values for temperature and gravity, we 
interpolated bilinearly in log T eff and log g within the Kurucz 
archive. 

Each pair of model and composite spectra was then con- 
volved with a Gaussian and the model regridded to the same 
wavelength scale as the corresponding composite. We then 
fitted the complete emergent spectrum from the model atmo- 
sphere to the composite, omitting only the region 7.43-11.25 



2278 COHEN ET AL.: SPECTRAL IRRADIANCE CALIBRATION 


2278 




Fig. 2. Eight 1.2-35 composites, smoothed by a Gaussian with FWHM of 0.3 /im, showing their ± 3 cr ranges (pairs of solid lines), compared with the 
identically smoothed relevant Kurucz models (short dashed lines), scaled optimally to match each composite spectrum, excluding only the region of the SiO 
fundamental (log X— 0.90). (These curves are to be read against the left-hand axes.) The major CO bands occur at log X— 0.67 (fundamental) and log X~0.36 
(first overtone) in these plots. The ratios of these smoothed spectra are also plotted in these figures (the dotted curves), in the form of composite divided by 
scaled model. These dotted ratio curves are associated with the right-hand axes. Note the SiO ratio absorption due to the absence of the lines of this 
fundamental band (log X—0.90) from the grid of models, and the apparent “emissions” in the CO fundamental (log X—0.67), where the models’ bands are 
too deep, (a) a Tau, (b) [3 Peg, (c) (i And, (d) f3 Gem, (e) a Hya, (f) a Boo, (g) y Dra, (h) a Cet. 
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Table 3. Fundamental parameters adopted for the generic Kurucz models.* 


Star 

Spectral Type 

T e // 

l <>g g 

a Tau 

K5III 

3898 

2.0 

P Peg 

M2.5II-III 

3600 

1-4 

P And 

MOIII 

3839 

1.5 

P Gem 

KOIII 

4844 

2.5 

a Hya 

K3II-III 

4141 

2.2 

a Boo 

K1III 

4362 

2.4 

7 Dra 

K5III 

3986 

2.0 

a Cet 

M1.5III 

3745 

1.3 


*This tale can be found in the AAS CD-ROM Series, Vol. 7, 1996. 


/xm, where the SiO fundamental is apparent in the observa- 
tions but lacking in the Kurucz models. The “fit” was ac- 
complished simply by treating the smoothed model and the 
smoothed composite as two independent spectral fragments 
that we wished to “splice.” This procedure (described in 
detail in Paper II) involves minimization of the sum- 
squared differences between the two spectral shapes, taking 
into account the variances of the observed points. The 
method provides the overall optimal scale factor for the 
model, and hence provides a value for the angular diameter. 
The results showed minimal sensitivity to sensible values of 
the FWHM used for the smoothing (i.e., in the range 0.1— 0.3 
fxm). 

Figure 2 illustrates the results of these direct comparisons, 
where we present the ratios of eight of our twelve composites 
(those with fullest wavelength coverage) to the optimally- 
scaled relevant Kurucz models. The fundamental parameters 
we used for the eight stars presented in these figures are 
summarized in Table 3. (The four spectra not shown present 
similar comparisons though the quality of fit is higher when 
we have a broader spectral range in the observations.) Apart 
from the SiO fundamental region (X—8.0 /zm), the models fit 
well, and the ratio spectra are essentially featureless (allow- 
ing for the 1 cr absolute uncertainties in our spectra), except 
for deviations of order 10% through the CO fundamental 
(\~4.7 /zm) and occasional, smaller, more ambiguous, fea- 
tures in the CO first overtone bands (\~2.3 /z m). We con- 
clude that the Kurucz models provide good predictors of 
stellar spectral energy distributions for these cool giants, at 
better than the ±2%-3% level outside the major CO and 
SiO absorption features, and about 5%— 10% in the CO 
bands. We recognize that a Boo is anomalous for its spectral 
type by virtue of its metal-poor character, yet this appears 
predominantly to affect the strong bands and not the less 
heavily blanketed regions where the spectrum approaches the 
true continuum. 

The discrepancies in CO might arise for a number of rea- 
sons, including our use of generic, as opposed to customized, 
models for these individual stars, Kurucz’ s use of compara- 
tively old CO line strengths (i.e., Kirby-Docken & Liu 
1978), and the omission of H 2 0 opacity in the Kurucz mod- 
els. Further investigation will be needed to sort this out. It is 
interesting to note that there is no clear tendency for stars of 
lower temperature to have systematically greater deviations 


Table 4. “Header” information accompanying new a Boo composite.* 

12-21-95 OBSERVED SPECTRUM OF ALPHA BOOTIS 
a Boo photometry file: photometry actually used to construct the spectrum 


Name 

FWHM 

(H 

Mag.iUnc. 

Eff Wvl 
(Vega) 
(pm) 

Eff Wvl 
(star) 
(fim) 

f a 

W cm -2 /zni -1 

Source 

Kn 

0.0488 

-3.07±0.01 

2.208 

2.205 

6.69E-13 

Selby, Hammersley 

Ln 

0.1443 

-3.15±0.01 

3.782 

3.762 

9.35E-14 

Selby, Hammersley 

M 

0.6677 

-2.97i0.02 

4.758 

4.738 

3.24E-14 

Strecker et al. 1979 

UK1RT87 

0.8611 

-3.13±0.01 

8.770 

8.779 

3.43E-15 

UKIRT Service data 

UKIRT98 

0.9455 

-3.16i0.01 

9.843 

9.828 

2.23E-15 

UKIRT Service data 

UKIRTll 

1.1198 

-3.16i0.01 

11.641 

11.639 

1.15E-15 

UKIRT Service data 

UKIRT12 

1.1782 

-3.19i0.01 

12.432 

12.427 

9.13E-16 

UKIRT Service data 


This UKIRT filter set represents a series of narrowband 10 fim passbands 
provided in common to UKIRT, Univ. Minnesota, and IRTF for ‘‘silicate" 
work. The “11" filter is centered near 11.7 n m, the“12" filter is centered 
near 12 . 5 fi m. 

Spectral fragments and portions of these actually used in observed spectrum 
(‘‘used" may include combination with other data where overlaps occur) 


Fragment 

Reference 

Total range 
(i‘ m ) 

Start and stop 
wavelengths (/tm) 

Average resolving 
power 

NIR 

1 

1.22-5.70 

1.22 5.58 

50 

KAO 

2 

3.65-9.39 

4.44-8.99 

150 

8-13 

3 

7.65-13.43 

7.65-13.43 

55 

LRS 

4 

7.80-22.70 

9.30-21.90 

30 

LONG 

5 

15.73 23.85 

15.73-23.85 

73 

VLONG 

6 

1.25-35.00 

21.80-35.00 

- 


References : 

• 1. Strecker, Erickson, ft Witteborn 1979, Ap.J. Suppl., 41, 501. 

• 2. FOGS data of May 11, 1992 KAO flight la Boo /a Lyr ] , and HIF0GS 
data of April 14, 1995 KAO flight [or Boo/c* Lyr] . 

• 3. FOGS Mt Lemmon data of Feb. 24, 1992 la Boo la Lyr], CGS3 UKIRT 
data of Hay 24 and 29, 1991 la Boo//? Peg], August 12, 1995 [a Boo/p 
Peg] and [a Boo / / 3 And] . The CGS3 spectra have the greatest weight 
in this combined data set. 

• 4. LRS raw data extracted from the new Groningen IRAS database and 
recalibrated with ‘‘LRSCAL" routine in ‘‘GIPSY" package. 

• 5. 20 fita UKIRT CGS3 data of May 24/25, 1991 for [a Boo//? Peg] . 

• 6. Engelke Fn. used for T=4362K (see Blackwell, Lynas-Gray, ft Petford 
1991, Aft A, 245, 567) and angular diameter=20 . 430 mas; we rescaled this 
to 21 . 12 mas . This Engelke Function was locked to the photometrically 
scaled combination of 8 — 13 and LRS spectra by splicing and used to 
replace the observations beyond 21.80 /zm. We applied an estimated 
uncertainty in EFn. of 2.4K, allowing for the change in shape of the 
EFn. for a temperature uncertainty of 100K at this effective temperature. 

INFORMATION ON SPLICES AND BIASES INCURRED 


Process 

Factor determined 

itBias (%) 

NIR cf. photometry 

1.003 

0.87 

813 cf. photometry 

1.031 

0.45 

LRS blue/ red bias 


0.03 

LRS splice to 813 

0.950 

0.02 

KAO joint splice to 
NIR and merged 813/LRS 

0.862 

0.41 

LONG splice to 
merged 813/LRS 

1.124 

0.66 

Engelke Fn. splice to 
combined 813/LRS /LONG 

1.017 

0.40 


*Actual spectra for other stars can be found in the AAS CD-ROM Series, 
Vol. 7, 1996. 

from expectation, in spite of the absence of H 2 0 lines from 
the Kurucz models. In the near future, we plan a much more 
extensive investigation of the influence of all the H 2 0, CO, 
and SiO lines that must be present in these giants on their 
emergent spectra. The overriding impression conveyed by 
Fig. 2, however, is of the strong predictive power of these 
models in representing these observed spectra. 

4. STELLAR ANGULAR DIAMETERS 

Table 5 (this table can also be found in the AAS CD- 
ROM Series, Vol. 7, 1996) summarizes four estimates of 
angular diameter for these stars: (i) from the multiplicative 
scale of the best- fitting Kurucz model appropriate to each 
star (as just described); (ii) from our Engelke approxima- 
tions; (iii) from the independent work of Blackwell et al. 
(1990), Blackwell et al. (1991) using the MARCS model at- 
mosphere code and a different absolute calibration for Vega 
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Table 5. Comparison of derived angular diameters with measured values.* 


Star 

Kurucz 

model 

Our Engelke 
function 

Blackwell 
et al. IRFM 

Observed 

Reference- 

or Tau 

20.92 

21.32 

20.63 

20.88 

Radgway et al. 1982 


±0.21 

±0.58 

±0.41 

±0.10 


/?Peg 

16.49 

16.98 

16.73 

16.75 

Di Benedetto & Rabbia 1987 


±0.16 

±0.51 

±0.33 

±0.24 


P And 

13.47 

13.71 

13.22 

13.81 

Mozurkewich et al. 1991 


±0.13 

±0.37 

±0.26 

±0.11 


(5 Gem 

8.01 

8.07 

7.87 

8.03 

Mozurkewich et al. 1991 


±0.08 

±0.17 

±0.16 

±0.08 


o Boo 

20.76 

21.12 

20.43 

20.95 

Di Benedetto & Foy 1986 


±0.21 

±0.51 

±0.41 

±0.20 


7 Dra 

10.04 

10.17 

10.00 

10.13 

Di Benedetto & Rabbia 1987 


±0.10 

±0.27 

±0.20 

±0.24 


o Cet 

12.52 

12.66 

12.64 

13.23 

Mozurkewich et al. 1991 


±0.13 

±0.36 

±0.25 

±0.20 


7 Cru 

25.49 

26.14 

- 

_ 



±0.25 

±0.86 

- 

- 


a Hya 

9.21 

9.27 

- 

- 



±0.09 

±0.24 

- 

- 


o TrA 

9.51 

9.55 

- 

- 



±0.10 

±0.26 


_ 


/i UMa 

8.16 

8.32 

- 

_ 



±0.08 

±0.28 

- 

- 


e. Car 

12.39 

12.67 

- 

- 



±0.12 

±0.33 

- 

- 



*This table can be found in the AAS CD-ROM Series, Vol. 7, 1996. 


than our own, or Blackwell & Lynas-Gray (1994), who re- 
vised their earlier results to use the Kurucz model grid and 
adopted our absolute calibration of the Selby narrowbands 
(Paper I); and (iv) from Michelson intensity interferometry. 
For the latter, we have used the ‘ ‘true’ ’ (limb-darkened) an- 
gular diameters (as opposed to those based on interpretation 
of the observed visibility fringes with the “uniform diam- 
eter” assumption). Details of the calculations used to convert 
observed to true angular diameters are presented by Di Bene- 
detto & Rabbia (1987), Mozurkewich et al. (1991), and in 
the more recent summary by Di Benedetto (1993). For a 
Tau, there is also an extensive series of measurements based 
upon lunar occultations (see White & Feierman 1987) with 
the best weighted value for this star being 20.88 mas (Ridg- 
way et al. 1982). (In the lower part of Table 5, we have 
included the five stars apparently without direct observa- 
tions, or IRFM-based estimates of diameter.) There is good 
accord between all four estimates for the seven stars for 
which all this information is available. The existence of such 
an accord makes a profound statement because each estimate 
is based on a very different spectral region. Our match to 
Kurucz models assesses the quality of fit over each entire 
composite (excluding the SiO fundamental), but this effec- 
tively assigns greatest weight to the shorter wavelengths (be- 
low ~7.5 /Am), where our total uncertainties are relatively 
small. By contrast, our Engelke approximations are fitted to 
the observed spectra longward of the SiO fundamental, typi- 
cally beyond 12 /Am. Blackwell and colleagues treat only 
narrowband photometry in the 1.2-3. 8 /Am region. Yet all 
these methodologies agree with the direct observations 
within the stated errors, and the technique of matching ge- 
neric Kurucz models to composites provides meaningful di- 
ameters and formal uncertainties in these diameters. 

The third column in Table 5 presents the most recent de- 
termination of angular diameter by Blackwell and colleagues 
for each of these stars. This corresponds to the Blackwell 
et al. (1991) paper for all stars except /3 Gem, for which 
Blackwell & Lynas-Gray (1994) provide a new value of 8.03 
mas with a methodology that is closest to our own (Kurucz 


models; our calibration of the Selby passbands). This value 
is essentially identical to our own diameter (8.01 mas). Un- 
fortunately, these authors have not published similarly reca- 
librated diameters for the other stars in Table 5. It is difficult 
to make accurate estimates of what values they would have 
obtained by application of their most recent method but, 
based on fl Gem, their diameters would have been quite 
close to our values. 

We note from Table 5 that the Kurucz models lead to 
systematically smaller angular diameters than indicated by 
the observations. We can speculate on possible causes of this 
phenomenon but we cannot fully account for it in the present 
paper. Several factors may contribute. First, the particular 
technique used here is designed to treat stars as though we 
lacked individual knowledge of their fundamental param- 
eters, including elemental and isotopic abundances (in readi- 
ness for application of the generic ‘ ‘template’ ’ scheme in the 
future). Perhaps, in this manner, or through our effective 
temperature- spectral class transformation, we have system- 
atically misrepresented the actual stars. Second, perhaps our 
attempt to fit the entire observable spectral region, without 
the SiO fundamental, has made our result more vulnerable to 
any inadequacies in the representation of the absorption 
strengths of the myriad molecular lines. (Recall in particular 
that the Kurucz and MARCS models do not include the po- 
tentially important H 2 0 opacity.) Use of a different scheme 
which weights the relative continuum portions of the spec- 
trum more heavily might lead to more secure diameters. 
However, even this is unclear since the fluxes in these ‘ ‘con- 
tinuum” regions are affected in reality both by a veil of 
molecular lines and by the brighter true continuum levels 
which arise when deeper atmospheric layers are heated by 
line blanketing. The Selby narrowband filters used to support 
the IRFM have isophotal wavelengths close to 1.25, 2.21, 
and 3.78 /im (logarithmic values of 0.10, 0.34, 0.58). These 
tend to correspond to local maxima close to unity in the ratio 
plots shown in Fig. 2. Therefore, basing our scaling of mod- 
els on the fluxes at these wavelengths might improve the 
agreement with the observations. 

The diameters implied by the Engelke approximations are 
generally bigger than those from using real models. This 
arises because Engelke functions are essentially Rayleigh— 
Jeans in character in the domain in which we employ them, 
whereas true model spectra decline more slowly with in- 
creasing wavelength. Therefore, to match the level of a given 
observed spectrum, one requires larger scale factors for En- 
gelke functions. 

It is worth emphasizing that the diameters derived from 
both the Kurucz and MARCS models generally accord, at 
the la level, with the observed diameters. This level of 
agreement implies discrepancies in diameter at only the 1% 
level, yet Blackwell & Lynas-Gray (1994) estimated conser- 
vative uncertainties in diameters derived by the IRFM of 4%. 
It is interesting that neither model grid reproduces the ob- 
served subtense of a Cet, but one must recall that angular 
diameters of resolved stars depend on limb darkening correc- 
tions to derive the “true” diameter from the observed. Such 
corrections might plausibly be suspect on very cool giants 
such as a Cet. An informative discussion of the complexity 
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Fig. 3. Ratio of calibrated models of Sirius and Vega (short dashed line), 
compared with CGS3 direct spectral observations (solid line), and with the 
measured broad and narrowband ratios (filled squares with ±lrr error bars) 
in well-characterized infrared photometric bandpasses. 

of limb-darkening in cool photospheres may be found in 
Scholtz & Takeda (1987). 

Finally, for y Cru, our angular diameter is quite different 
from the 41.1 mas cited by Judge & Stencel (1991). How- 
ever, their value was an. indirect one, based upon adopting 
colors similar to those of other very cool M giants and ap- 
plying relations due to Barnes et al. (1978) between color 
and visual surface flux (Stencel 1996; cf. Judge & Stencel 
1991, Appendix A, Sec. 2, last sentence). 

5. THE SPECTRAL RATIO OF SIRIUS/VEGA 

The derivation of plausible angular diameters from our 
composites provides external validation of their accuracies. 
This, in turn, vindicates the procedures we have devised to 
assemble composite spectra. The derivation (Paper VI) of an 
angular diameter for a 1 Cen, that is consistent with the lit- 
erature for that star, is a further indirect check on the 
adopted, calibrated model for Sirius. As a final cross check 
on the overall scheme we have sought data not used in Paper 
I on the ratio of Sirius to Vega. Two kinds of data are of 
value: independent photometric measurements of Sirius 
compared with Vega; and the ratio of their spectra. 

Our calibration of the Kurucz Sirius model relative to the 
calibrated model for Vega (Paper I) was based upon the bo- 
lometry of Cohen & Barlow, published by Deacon (1991). 
Therefore, we sought other well-characterized photometric 
systems for which the difference in magnitude between these 
hot stars has been measured. We require detailed cold- 
scanned filter transmission profiles for these systems, and we 
have applied an atmosphere appropriate to each site (cf. Pa- 
per I). Of particular value are data by Carter (1996: detailed 
in Paper VIII) in the SAAO JHKL system [as characterized 
for us by Glass (private communication)], and by Hammer- 
sley (1995: JHK; to be represented in detail in Paper VIII) 
and Alonso et al. (1994: JHKL') in the Tenerife broadband 
InSb passbands. The JHK data of Alonso et al. are already 
included in Hammersley’s reanalysis of the entire photomet- 



Fig. 4. As for Fig. 3, but now combining direct and indirect CGS3 spec- 
troscopy of Sirius/Vega, and illustrating the ±2 cr bounds on the ratio of the 
model spectra. 

ric archive of the Telescopio Carlos Sanchez (Paper VIII) so 
only these authors’ L' result is separately represented here. 
Additional photometry can be found in S inton & Tittemore 
(1984) using L' and M passbands which, from their descrip- 
tions, are identical with the UKIRT bolometer passbands of 
the same name that we calibrated in Paper I (although Sinton 
and Tittemore’ s measurements were never used in Paper I). 
We emphasize that we have taken great pains to avoid cir- 
cularity: we do not use our expectation for Sirius’ magni- 
tudes in these passbands from a recalibration of the authors’ 
systems, but focus strictly on the observed differences in 
magnitude as measured by these authors. Figure 3 illustrates 
the ratios equivalent to each of these magnitude differences 
between Sirius and Vega. Filled squares represent the pho- 
tometry with ±lcr error bars, defined as the root-sum-square 
of the observed magnitude uncertainties of these two stars, in 
each passband. Points are plotted at their isophotal wave- 
lengths (almost identical in these bands for Sirius and Vega: 
cf. Paper I). The longest wavelength point plotted in Fig. 3 
and 4 is the high precision measurement made by IRAS at 12 
/xm (from Table VI.C.3 in the IRAS Explanatory Supplement 
1988). 

We wish to make a comparison of our adopted ratio of the 
two model spectra, and to assign meaningful uncertainties to 
their ratio. The relevant uncertainty for the spectral ratio of 
Sirius to Vega is the uncorrelated error, rather than the un- 
derpinning 1.45% absolute uncertainty in Vega’s calibration 
which is directly transferred, as a bias (a correlated uncer- 
tainty), when we calibrated Sirius relative to Vega. If we 
assign uncertainties of ±0.01 111 to each of Deacon’s magni- 
tudes, then the root-sum-square of these random errors with 
the small bias associated with locking the Sirius model to the 
expected photometry of Deacon, yields a conservative esti- 
mate of ±0.5%. Figure 4 compares this photometry with the 
±2cr bounds in the mean ratio of our adopted calibrated 
spectra for these stars (short dashed lines). The photometry 
points appear to follow the upward trend (with decreasing 
wavelength) in the models’ spectral ratio very well. This is 
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Fig. 5. Further evidence for the Rayleigh -Jeans character of cool stellar 
spectra in the form of spectral “flatness” in (\,\ 4 F\) space. UKIRT CGS3 
16-24 /an spectra for a Boo (unsealed, solid line) and a Cet (scaled by a 
factor 3.2, dash-dot line) are compared with their relevant Engelke approxi- 
mations (dashed lines across the middle of the plot). That for a Boo appears 
slightly above that for a Cet. 

entirely due to differences in the two adopted stellar effective 
temperatures. 

What observational spectroscopic evidence can we ad- 
duce in support of our models’ ratio? There are two kinds of 
measurement to discuss, because it is a challenging observa- 
tion to secure spectra of Sirius and Vega in a single night. 
Given data on the two stars in one night, one can construct 
the direct spectral ratio. The second, indirect, measurement 
of spectral ratio that we consider is formed by use of inter- 
mediate stars, so that one compares an airmass-matched ratio 
(cf. Papers II, V) of Sirius to a bright star, and of Vega to the 
same bright star, also matched in the airmasses of observa- 
tion. The indirect method offers the advantage of being able 
to choose good bright intermediaries and to secure close 
matches in airmass, at the cost of added noise in the eventual 
ratio of Sirius to Vega, attendant on the products or quotients 
of these intermediate ratio spectra. 

We have found one direct measurement of Sirius to Vega 
from the UKIRT archives. This is constructed from 10 /zm 
spectra taken with the CGS3 spectrometer, by Skinner and 
Sylvester on 1992 October 4, and kindly furnished for our 
use by Sylvester (1996). The two stars were not closely 
matched in airmass (Vega z = 1.29; Sirius z = 1.55). Conse- 
quently, the short-wavelength edge was lost to water vapor 
differences and the 9.6 /xm telluric ozone feature was uncan- 
celed. This necessitated deletion of the affected data (7.5-7. 7 
and 9.75-9.9 gm). The resulting spectrum was normalized 
using photometry from Deacon (1991) in the UKIRT 8.7 yum 
“narrowband” filter. Figure 3 presents this direct spectral 
ratio, as the simplest comparison we can make with the mod- 
els’ ratio. 

To enhance signal to noise, we needed to secure addi- 
tional indirect spectral ratios, so we sought airmass matching 
with bright intermediaries. Two such indirect routes are 
available to us from the UKIRT CGS3 archives. One in- 
volves the quotient of ratios of spectra, [/3 Peg /a Lyr] (1995 


Fig. 6. Far-infrared model continuum extrapolation (solid line) of our a Tau 
composite, compared with the extrapolated Engelke approximation (short 
dashed line) for this star. The maximum difference between these two spec- 
tra occurs near 65 gm, where it attains almost 2%. 

August 12) and [/? Peg la Tau] (1993 November 4 and 5), 
from which we derive an estimate of the spectral ratio [a 
Tau/tr Lyr], The second similarly involves [/3 And/a; Lyr] 
(1995 August 11 and 12) and [y8 And/or Tau] (1993 August 
28), whose quotient gives another estimate of [a Tau / a Lyr]. 
We combined these two ratio spectra (with inverse-variance 
weighting) to yield a cleaner [a Tau / a Lyr], which we di- 
vided by our published ratio of [or Tau /a CMa] (from 1991 
November 9: Fig. 4 of Paper III), to create [or CMa/ or Lyr]. 
Finally, we combined this indirect ratio (which has broader 
wavelength coverage than the direct method because of the 
close matching in airmasses of the quotient spectra) with the 
direct approach, using inverse variance weighting. Essen- 
tially the same spectral ratio results as from the direct mea- 
surement alone, but now it is seen with higher signal to noise 
and over a larger and more completely sampled wavelength 
interval. The resulting ratio spectrum appears in Fig. 4, 
where we have indicated the ±2<x bounds on this observed 
spectrum, exactly as shown for the ratio of model spectra. 
The two ±2cr swathes (model and CGS3) overlap signifi- 
cantly. 

Therefore, we conclude that both independent photometry 
and spectroscopy support the ratio of our two, independent, 
adopted models for the spectra of Sirius and Vega. 

6. THE BEHAVIOR OF REAL GIANT STARS FROM 10 TO 35 /xm 

The agreement between cols. (1) and (2) of Table 5 sug- 
gests that the simplistic Engelke function is a reasonable 
approximation to the general shape of emergent spectra pre- 
dicted by fully detailed model atmosphere calculations. 
Closer inspection reveals that the diameter implied by the 
Engelke approximations is invariably larger than that from 
matching to real models. This happens because Engelke 
functions decline slightly faster with wavelength than model 
spectra, and therefore require slightly greater multipliers to 
match a given observed spectrum. What of real stars? 
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Only two normal giant stars, a Tau and y Cru, were ob- 
served by Glaccum (1990) from the KAO , and these spectral 
fragments owe their shapes to standard planetary absolute 
calibration techniques (Paper II; Moseley et al. 1989). These 
two stars do seem to behave in a Rayleigh-Jeans fashion 
once the spectrum rises out of the strong SiO fundamental 
absorption; i.e., in \ 4 F\ space, their observed spectra are 
essentially independent of wavelength beyond about 12 yum. 
Additional support for the observed “flatness” of real stellar 
spectra in this space (Figs. 1 and 2) comes from our CGS3 
20 /xm spectra taken at the UKIRT and, to a lesser extent 
(due to the noise at the longest wavelengths), from the LRS 
too (cf. for a Tau in Fig. 11 and 12 of Papers II; for y Cru, 
Fig. 5, and (3 And, Fig. 6, both in Paper IV). In the present 
paper, we adduce additional evidence for the Rayleigh-Jeans 
nature of real stellar spectra from 20 /xm CGS3 spectra of a 
Boo and a Cet. These are compared, in Fig. 5, with the 
relevant Engelke approximations. The spectral “flatness” is 
emphasized by these direct plots of the calibrated stellar 
spectra in X 4 F k space. 

7. SYNTHESIZING IRAS 12/25 FLUX DENSITIES 

We hoped further to probe this apparent and persistent 
discrepancy between the long-wavelength (beyond —12 /xm) 
behavior of real and model stars, by comparing our compos- 
ite spectra with IRAS observations. Every composite is fun- 
damentally built from observations based on our model spec- 
tra for Sirius and Vega. Therefore, we would expect to 
conclude in favor of the same “recalibration” factors for 
IRAS flux densities deduced from photometry alone in Paper 
I, by use of IRAS flux densities synthesized from our com- 
posites, unless real stars deviated significantly from 
Rayleigh-Jeans slopes. The procedure is straightforward, re- 
quiring only the integration of the IRAS system response 
functions for the 12 and 25 /xm passbands over our spectral 
composites, and direct comparison with the IRAS in-band 
fluxes, and the flux densities derived from these upon divi- 
sion by the bandwidths (cf. Paper I). 

Table 6 (also found in the AAS-CD-ROM Series, Vol. 7, 
1996) summarizes the results of these direct comparisons 
with the IRAS Faint Source Survey (Moshir et al. 1992), re- 
putedly of higher quality than the IRAS Point Source Catalog 
(1988: hereafter referred to as PSC) fluxes, and with the PSC 
itself. Table 6 incorporates both the uncertainties that attend 
our derived fluxes and those inherent in the IRAS measure- 
ments (“RELUNC” quantities in these catalogs). The com- 
parison suggests that these IRAS point source products over- 
estimate 12 and 25 /xm flux densities by statistically the 
same amount in the PSC and the FSS, namely, by 4.1% at 12 

Table 6. Ratios of IRAS flux densities synthesized by integration of the 
IRAS passbands over our observed composite spectra (12,25 /im), extrapo- 
lated by model atmospheres (60,100 /tm), to those actually observed by 
IRAS.* 


IRAS data 

12 urn 

(T 

25 

a 

60 //,m 

<7 

100 fim 

<T 

FSS 

0.961 

0.012 

0.946 

0.013 

0.997 

0.015 

1.038 

0.040 

PSC 

0.961 

0.010 

0.943 

0.010 

0.985 

0.029 

1.015 

0.026 


This table can be found in the AAS CD-ROM Series, Vol. 7, 1996. 


and by 5.7% at 25 /xm. From this we conclude that our 
spectral composites, based on the calibrated spectra of Sirius 
and Vega (Paper I), represent self-consistent products be- 
cause these factors for IRAS’ s calibration, derived here from 
cool stellar spectra (0.96±0.01 and 0.94±0.01), are in ac- 
cord (at the 1 cr level) with those determined purely from our 
basis models for Sirius and Vega (12 /xm: 0.98±0.01; 25: 
0.94±0.02). It is also tempting to conclude that, in the en- 
semble average, IRAS measurements vindicate the essentially 
Rayleigh-Jeans character of the infrared photospheric emis- 
sion from these cool giants in the 7-35 /xm range, although 
the actual IRAS uncertainties on any star are too large to 
reach this conclusion for an individual object. 


8. EXTRAPOLATION OF COMPOSITES TO 300 /xm 

The 35 /xm limit of our standard composite spectra does 
not afford the same test with longer wavelength IRAS data. 
Consequently, we have adopted the emergent spectral shapes 
of a grid of mid- to far-infrared continuum models of red 
giant branch stars, computed as described below. We inter- 
polated in this grid to provide a model that matched the 
correct effective temperature for each composite (and we as- 
sumed the red giant branch gravities interpolated were cor- 
rect for our observed giants). The attributes of the giants 
come from VandenBerg & Laskarides (1987) and are close 
to the expected gravities for our giants. By fitting these 
shapes to the observed portions of composites, longward of 
the SiO fundamental’s influence, but shortward of the wave- 
length where we chose to replace noisy data by Engelke 
functions, we created best-fitting extensions of each compos- 
ite. 

In order to obtain representative effective temperature- 
log(gravity) (T efr ,log g) combinations for the late-K and 
early-M giant branch stars of the disk in our observed 
sample, we adopted the theoretical red giant branch tracks 
computed by VandenBerg & Laskarides (1987). We chose 
the particular track with the following model parameters: 
^/^©=1.3; Y =0.25; Z-0.0169; [M/H]=0.0. 

Stellar evolution models with these parameters have an 
age of 5.28 Gyr at the He flash and, thus, should be reason- 
ably representative of evolved stars of the disk. To obtain the 
(T eff ,log g ) needed for the stellar model atmosphere calcula- 
tions, we selected three combinations of (T eff Jog g), spaced 
by 500 K, on the aforementioned VandenBerg & Laskarides 
track: (4500,2.86), (4000,1.87), and (3500,0.975). 

To obtain a model which might be appropriate for the 
hottest star in our sample, (3 Gem, we adopted the specific 
parameters of T eff =4844 K and log g — 2.24, which are in 
excellent agreement with the analyses of this star by Black- 
well & Lynas-Gray (1994), Di Benedetto (1993), and Ruland 
et al (1980). 

To obtain the temperature vs. continuum optical depth at 
1 /xm relations [T , (tj m )], needed to compute emergent 
fluxes for our four desired models, we began with published 
model grids and proceeded in the following fashion for each 
of the desired models: 

(1) the temperature-pressure relation of the grid model 
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nearest in 7'eff and log g to the desired model was used to 
obtain the grid model’s T(t\ ); 

(2) the grid model’s T(t c { Mm ) relation was then scaled by 
the ratio 7’ eff (desired)/r eff (grid model). The validity of this 
procedure over moderate ranges in T eff has been demon- 
strated by Carbon & Gingerich (1969) and Gustafsson et al. 
(1975); 

(3) and, finally, using the resultant T{r c x /U>m ), the equation 
of hydrostatic equilibrium was reintegrated using the log g of 
the desired model. 

All the calculations were carried out using the set of solar 
abundances adopted by the authors of the respective model 
grids. 

For the three coolest models, (4500,2.86), (400,1.87), and 
(3500,0.975), we interpolated in the model atmosphere grid 
published by Brown et al. (1989). For the yS Gem model, 
(4844,2.24), we interpolated in Kurucz’s (1991, 1993) model 
grid to obtain the required T— P relation. For each of the four 
models, we computed continuum spectral energy distribu- 
tions at 101 wavelengths from 1.0 to 300 /x m using the 
SOURCE model atmosphere program described in Carbon 
et al. (1982). For ISO purposes, we created products gridded 
every 5 /xm from 25 to 300 /xm to supplement and extend 
our absolute calibrated composites out to the longest wave- 
lengths relevant to ISO. 

We assessed uncertainties in these extrapolations based 
on the recent study of the sensitivity of far-infrared continua 
to fundamental stellar parameters by van der Blieck et al. 
(1996). Explicitly, at 100 /xm, these are ~ 2 % (corresponding 
to ±100 K in effective temperature), ~0.2% (±0.3 in log g ), 
and *£l% (±0.2 in metallicity). In fact, as these authors 
point out, errors in gravity are even greater for late-type gi- 
ants than for ensembles of warmer stars, typically attaining 
±0.5 in log g. Van Blieck et al. (1996) also simulated chro- 
mospheres and considered the influence of circumstellar 
dust. We examined their estimates for uncertainties in far- 
infrared flux due to errors in fundamental parameters, in tem- 
perature structure, and other causes, and adopted 4.8% as 
typical for the root-sum-square combination of all these 
components for the mid-G to early-M giants they considered. 
To these we have appended an extra 3% to account for the 
typical expected contributions to the spectra of molecular 
line blanketing due to the copious lines of water vapor, CO, 
and SiO, and their isotopes. 

9. SYNTHESIZING IRAS 60/100 /mi FLUX DENSITIES 

These extrapolated composites permit integration over the 
60 and 100 /xm IRAS bands. Table 6 presents the results of 
direct comparison of long- wavelength flux densities synthe- 
sized from our spectra with those measured by IRAS. In spite 
of the rather larger relative uncertainties, “RELUNCs,” as- 
sociated with the IRAS long-wavelength measurements we 
can, nevertheless, determine from this ensemble of compos- 
ites that IRAS measurements are consistent with our cali- 
brated observed spectra, extrapolated through theoretical 
model atmospheres. 

It is, perhaps, also of interest to examine the long- 
wavelength quality of the simplistic Engelke approximations. 


The absolute difference in emergent spectra of these giants, 
between their best-fit Engelke approximations and our grid 
of far-infrared continua, is always small. In fact, as Fig. 6 
illustrates, for a Tau it amounts to a mere 2% over the range 
25-200 /xm. For the set of 12 composites now available, this 
maximum difference between Engelke approximations and 
model continua lies in the range 2%-5%, and always occurs 
between 60 and 70 /xm. As a further test, we also contin- 
uedthe same Engelke functions that best fitted each compos- 
ite (see Papers IV, VI, and this paper) to wavelengths long 
enough to synthesize IRAS flux densities by integrating these 
extended spectra over the long IRAS system passbands. The 
corresponding ratios of synthetic 60 and 100 /xm flux densi- 
ties predicted by Engelke functions to those observed by 
IRAS (the FSS) are: 0.964±0.013 (60 /xm) and 0.976±0.037 
(100 /xm). These should be compared with the corresponding 
numbers in the first line of Table 6. Differences are only at 
the 1 .6er level, so we could not discriminate statistically be- 
tween Engelke approximations and true models of the infra- 
red continuum solely on the basis of IRAS flux densities. The 
fundamental limitation again lies in the “RELUNCs.” 

10. CONCLUSIONS 

We have assembled absolutely calibrated, complete, con- 
tinuous stellar spectra for five cool stars: a Boo (updated), y 
Dra, a Cet, y Cru, and /x UMa. These spectra are constrained 
by carefully characterized photometry and have an estimated 
absolute uncertainty (within our adopted absolute frame- 
work: Paper I) of order 3% (la) across most of the 1.2-35 
/xm range. Such spectra establish the pedigree for secondary 
stellar standards and thereby the flexibility to calibrate arbi- 
trary filter systems and sensors on the ground, from air- 
planes, or satellites. We find good agreement in shape be- 
tween these spectra and generic Kurucz model atmosphere 
spectra, with the obvious exception of the SiO fundamental 
and some systematic small (~10%) deviations associated 
with the CO fundamental. Stellar angular diameters derived 
from our spectra are in accord with directly measured diam- 
eters. 

Using a grid of computed stellar model continua, we have 
extrapolated our observed spectra as far as 300 /xm, main- 
taining our absolute context, to provide calibration support 
for the ISO PHT and LWS instruments. Tests of these ex- 
trapolations using IRAS 60 and 100 /xm point source flux 
densities show consistency with our predictions. The com- 
posite spectra with their associated pedigree files will appear 
on an AIP Journal CD-ROM and eventually will be obtain- 
able through NASA’s NSSDC at Goddard Space Flight Cen- 
ter. 
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20) jjm CGS3 spectra that constitute a vital part of our ob- 
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We thank Dr. P. Wesselius of SRON, Groningen, for provid- 
ing NAS A- Ames with its copy of the “LRSVAX” archive 
which constitutes a vital resource in our ongoing effort in 
infrared calibration. 
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